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A. BACKGROUND 


For more than two decades, the need to solve control 
problems in aerospace application has been the primary 
driving force behind the modern control theory development. 
Problems in manoeuvring, guidance and tracking of aircrafts 
and space vehicles have motivated the development of various 
control design and synthesis methods. One of the methods, 
the so-called Linear Quadratic Control or LQ-control 
[Refs. 1,2] has been widely used and is treated extensively 
in the control literatures. Unlike most classical methods 
where the design are based on conventional time response 
criteria, the LQ theory treats the problem of designing 
memerollers as that of minimizing a quadratic cost function 
of states and control inputs. The design problems become 
that of selecting suitable weighting matrices in the 
performance index. Two questions naturally arise from this 
method. First, how does one select the weighting matrices 
and second, once a set weighting matrices is selected, how 
does one know that it is a good design. There are generally 
two approaches to the first problem; the obvious one is to 


rely on physical arguments and a certain amount of trial and 


EN or (Refs. 3,4] Unfortunately, such formulation can be 
obtained in only a few cases. Reference 4 provides a few 
guidelines that can be employed. The second approach is to 


avoid the physical aspect of the performance index but 
instead try to relate the weighting matrices with some other 
performance specification. For example, Tyler and Tuteur 
[Ref. 5] expressed the characteristic polynomial as ап 


explicit function of the weighting matrices for single input 


1 


single output (SISO) system with diagonal weighting 
elements. Root-Locus type of procedure were used to show 
how the variation in weighting elements affect the eigenva- 
lues of the closed-loop system. Similar relations were 
explored in [Refs. 6,7] in which more general expressions 
were obtained. Their uses were restricted to single input 
case due to difficulty in handling polynomial matrices. 
Solheim [Ref. 8] later developed a sequential design proce- 
dure based on diagonalized (decoupled) system. More recent 
results of eigenvalue placement in optimal control problem 
are presented in [Refs. 9,10,11,12]. 

The second issue of LQ design, 1.e., whether a good 
design has been obtained once the performance index has been 
fixed, is related to the multivariable nature of system. 
Unlike the single input case where the closed-loop eigenva- 
lues uniquely define the feedback gain and hence the 
weighting matrices, the MIMO structure provides additional 
dimension which allows further tradeoff for properties other 
than the closed-loop eigenvalue location. An example is the 
gain and phase margin, or in MINO Ease, the so-called 
robustness criteria. Robust-control has been the subject of 
extensive researches in recent years and results relating to 
optimal control can be found in (Refs. MEI 

It becomes evident from the above discussion that LQ 
control design and synthesis are not just a matter of speci- 
fying performance index. Optimal in the sense of satisfying 
performance index and perhaps pole location does not neces- 
sary means that a good design has been obtained. Other 
criteria like disturbance rejection, robustness and sensi- 
tivity need to be considered and incorporated in the design 
procedure. This is the motivation behind our present 
research which in turn leads to the development of the 
synthesis package. An overview of the thesis is described 


in the next sees tone 


[2 


B. OVERVIEW 


Im this section an overview of the thesis is given. А 
background of multivariable optimal control theory in terms 
of its structure, frequency domain characteristic and asymp- 
totic properties is first discussed together with derivation 
of some useful relationships in Chapter 2. General robust- 
ness concepts and its application in Linear Quadratic (LQ) 
Control is presented in Chapter 3. In Chapter 4, a computer 
aided design package for pole-placement synthesis based on 
numerical optimization technique is presented. This package 
provides a useful computational tool to support the material 
in the remaining chapter. A step-by-step pole-placement 
Synthesis procedure is also presented to illustrate how the 
package can be used to design optimal control system that 
meet time response criteria and other properties. The use of 
the pole-placement synthesis package in two actual design 
problems is demonstrated in Chapter 5. Results in terms of 
frequency and time response properties ,robustness (singular 
value decomposition) are compared to those obtained by other 
design methods. Program listings and an example of a design 


session are included in the Appendix. 
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II. MULTIVARIABLE LINEAR QUADRATIC CONTROL 


In this Chapter a brief review of Multivariable Linear 
Quadratic Control theory 45 presented. Structure of LQ 
system is first given, general stability properties are then 
discussed. In section two the general pole assignment 
problem is formulated for the MIMO state feedback system. 
Interpretation of the state and control input weighting 
matrices and their effect on the closed-loop time response 


behaviors and the asymptotic properties are presented. 


A. LINEAR QUADRATIC SYSTEM 


Consider the following linear time-invariant state-space 


System given by the equations: 


x(t) Ax(t) * Bu(t) (eqn 2.1) 


y(t) Cx(t) (eqn 2728 


where x(t) is the n-dimensional state vector, u(t) is the 


m-dimensional control vector and y(t) is the p-dimensional 


output vector. A ,B and C are real constant matrices ٣ 
dimension nxn, nxm and pxn respectively. Assuming that they - 
form a controllable pair(A,B) and an observable pair(A,;€)8 


the optimal feedback control law is obtained by minimizing 


the following quadratic performance index. 


Ca 
tt 


со 
j ( x! Q x + pulR u)dt (eqn 2.3) 
0 


where R is positive definite (R>0) for bounded input and p 


is a scalam Q is a semi-positive definite matrix (Q20). 
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When both Q and R are diagonal matrices, p defines the 
relative weight between the state and control weighting 
matrices in the performance index. The quadratic form x'Qx 
and u'Ru provide a weighted measure of the magnitude of the 
states and control vector respectively. The steady state 


control law that minimizes J is given by, 


mre) = -Fx(t) (eqn 2.4) 


where F is the feedback gain matrix which is given by 
=-=} T 
F= -RBP (eqn 2.5) 
The positive definite matrix P is given by the solution to 


the steady state Riccati equation. 


PC + AP + О- PBR'B'P = 0 (eqn 2.6) 


Equation 2.5 and 2.6 are well-known results in optimal 


control theory that yield the optimal closed-loop system 


x(t) = [A + BF]x(t) (eqn 2.7) 


whose closed-loop poles are given by the eigenvalues of the 
matrix [A+BF]. The LQ system given by equations 2.1 through 
2.6 is closed-loop stable and сап Бе represented by the 


feedback configuration in both the time and frequency domain 


as shown in Figure 2.1 
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Frequency Domain 


Time Domain 


Figure 2.1 State Feedback System-Time and Frequency Domain. 
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In Figure 2.1, if the loop is broken at the input (X) as 


shown, the loop transfer function is given by 


G(s) = F(sI - A) B (eqn 2.8) 


The matrix [I + G(s)] is called the Return Difference Matrix 
and will be shown to have some important feedback properties 


in the following section. 


B. POLE ASSIGNMENT PROBLEM 


In the most general form, the state feedback pole 
assignment problem in control system design can be formu- 


lated precisely as follows: 


"Given real matrices (A,B) of order (nxn, nxm) respec- 
tively and a set of n complex numbers (23,,24,,...44,) closed 
under complex conjugation. Find a real mxn matrix F such 
anaie the eigenvalues of [A+BF] areà;, i= 1, 2,3... n." 


In general, the closed-loop eigenvalues of (A+BF) can be 


arbitrary located in the complex plane, with the only 
restriction that complex characteristic eigenvalues must 
occur in complex conjugate pair. In other words, if the 
Matrices (A,B) are a completely controllable pair, the 


stability of the system can always be improved by state 
feedback. If the (A,B) pair is not completely controllable, 
then the system is required to be 'stabilizeable' meaning 
that in its controlability canonical form given by equation 
2.9 below, the matrix Any is asymptotically stable or any 
unstable subspace of equation 209 Xs also in the 


controllable subspace. 
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x, (t) P А , xa b B аа 
کت‎ О ТА >x, (Om o ٣ 
(eqn 2.99 


The solution to the pole assignment problem in the 


single input (m-1) case, when it exits, can be shown to be 


unique. In the multiple input case (l<m<n), the solution of 
the so-called inverse eigenvalue problem is, in general, 
undetermined with many degrees of freedom. Additional 


conditions must be Supplied in order to eliminate the extra 
degrees of freedom. This has been an area of active 
research in recent years anda number of approaches have 
been developed to relate the extra degrees of freedom with 
properties such as system elgenvectors, transmission zeros 
and robustness. (Reese. don) ie In this work, it is ۶٦ 
that Linear Quadratic formulation incorporating equation m EE 
will partly remove the uncertainty that exists in the 
multiple inputs case. It will be shown that useful proper- 
ties like robustness are guaranteed. It is also shown that 
the LQ type of pole placement formulation, when combined 
with eigenvectors assignment type of formulation, will 
produce some very useful design and synthesis procedures. 
In the present work, however, only the LQ eigenvalue place- 
ment problem is addressed; the problem of combined eigen- 


value and eigenvector assignment is briefly described in 
Chapter 5. 
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C. WEIGHTING MATRICES AND SYSTEM PERFORMANCES 


This section briefly reviews the effect of weighting 
matrices on system performance. The physical aspects of the 
weighting matrix for both single and multiple input systems 
are presented first together with some discussion on the 
asymptotic behavior of the LQ system. 

It was shown in the last section that under complete 
controllable conditions, the time-invariant linear system 
can be stabilized by a linear feedback control law. For the 
regulator-type problem where the aim is to bring the system 
from an arbitrary initial state to the zero state, the 
closed loop poles can be chosen far to the left on the 
complex plane. Convergence to the Zero state is fast but 
the large input required may not be practical. This natu- 
rally leads to an optimization problem where the trade off 
15 between speed of convergence to zero and the magnitude of 
the input. This is reflected in the two quadratic terms in 
the performance index. The quantity x'Qx in the first term 
of the performance index is a measure of the extent in which 
the state at time t deviates from the zero state. The 
matrix Q determines how much weight is attached to each of 
the component of the state. The integral J(xTQx)dt 15 а 
criteria for the cumulative deviation of x(t) from the zero 
State during the interval. 

The problem of large control input is resolved by 
incooperating the second quadratic term, J(u' Ru)dt. Larger 
value in the element of the control weighting matrix R will 
But in smaller input. It will be shown later that one can 
manipulate R to achieve some secondary design objectives. 
The remaining parameter وم‎ that need to be specified, 
accounts for the relative weighting between the state and 
the control inputs. Selecting optimal value for p depends 


on the particular problem and the design requirement. AS p 


1 


decreases, the integrated square regulating error decreases 
but the integrated square input increases. Very often, the 
optimal control problem is solved for many different values 
of p `, A graph similar to that shown in figure M 
obtained where the integrated square regulating error is 
plotted versus the integrated square input. An appropriate 
value of p is then selected to give sufficiently small 


regulator error without excessively large input. 


e increasing 
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Figure 2.2 Selecting Relative Weighting Parameter p. 
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The special case where p approaches zero ,the so-called 
asymptotic properties, has been shown to provide good 
insight for LQ control system design. Some of the results 
from [Refs. 4,18] are summarized below. 

As p decreases to zero for the system given by equa- 
tions 2.1 through 2.6, some ( say q ) of the closed-loop 
poles go to infinity while other ( n-q ) stay finite. Those 
remaining finite approach the left half plane zero of 
det[B'(-sr-A' )^'Q (sI A) 7 If m is the dimension of the 
input vector, then at least m closed-loop poles approaches 
1٦1 716127۰ A11 closed-loop poles that go to infinity do so 
by grouping into several Butterworth patterns of different 
order and radii. For p approaching infinity, the closed- 
loop poles approach the mirror image of the plant open-loop 
poles. To illustrate the asymptotic concept, two examples 
from [Ref. 4] are given below. 

An example of the asymptotic properties of a single 
input system is shown in Figure 2.3 . The closed-loop poles 
of a position control system using LQ type feedback are 
plotted as a function of P . The system has two open-loop 
poles at -4.6 and 0.0 . As p decreases, the closed-loop 
ئ٢‎ go to infinity along two straight lines that make an 
angle of T/4 with the negative real axis. AS p 
approaches infinity, both closed-loop poles first meet on 
the negative real axis and then approach the’ open-loop 
poles at (-4.6, 0). 

Figure 2.4 shows the asymptotic loci of the closed-loop 
poles for a multiple input system where n=4 and m=2. There 
are four open-loop poles at (-0.006123, +j0.09353) and 
120250. —*31.394) and both the state and control weighting 


matrices are taken to be of a diagonal form. AS p 
approaches to zero, опе closed-loop pole stay finite at 
-1.002. The remaining three go to infinity, two of which 


assume a second order Butterworth pattern and the last one 
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approaches on the negative real axis. When p approaches 
700 ИН all closed-loop poles approach the open-loop 
poles. 

Many researchers have explored this asymptotic proper- 
ties. Stein [Refs. 9,10] developed a procedure to select 
the control weighting matrices for a desired closed-loop 
asymptotic eigenstructure. It will be shown in the subseq- 
uent Chapters that the asymptotic properties provide useful 


guideline for the design procedures to be described. 


p -0 00002 





Figure 10 Single Input Asymptotic Root Loci. 


Д2 





Faraway Poles 





Nearby Poles 


Figure 2.4 Multiple Input Asymptotic Root Loci. 
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III. ROBUSTNESS THEORY AND LINEAR QUADRATIC CONTROL 


In the last Chapter, the effect of weighting matrices mam 
the closed-loop pole of the LQ system was discussed. This 
Chapter addresses yet another important feedback property 
that Col rol system designers are concerned wit 
Robustness. Robustness theory was developed when it was 
realized that the classical single loop Nyquist test was not 
adequate to guarantee stability when the multivariable open 
loop plant deviates from its model due to a variety of 
reasons [Ref. 15]. In the following section , the concept 
of robust design for both SISO and MIMO general feedback 
System 1S reviewed. The Singular value analysis is 
discussed in term of multiplicative type of disturbances. 
Finally, robustness for linear quadratic state feedback are 
presented in terms of the effect of weighting matrices on 


the singular value curve. 


A. ROBUSTNESS CONCEPTS 


The concept of robustness for SISO system can best ٣7 
described in terms of the definition of phase and gain 
margin. A system characterized by good gain and phase 
margin implies that changes in the plant model parameters 
and changes in the loop gain and/or phase may be accommo- 
dated without loss of stability. The gain and phase margins 
of a SISO system are defined with reference to the perturbed 
System in Figure 3.1 

Assuming that the unperturbed system (l(jw)=1) is 
stable, the positive (or negative) phase margin is the value 
of ẹ greater (or less) than zero at which the perturbed 


system with l(jw) - k*exp(je) becomes unstable. The upward 
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Figure 3.1 The Perturbed SISO System. 


(or downward) gain margin is the smallest (or greatest) 
value of l(jw)=k=constant, for which the system become 
unstable. With reference to the classical Nyquist plot in 
Figure 3.2 , gain and phase margin are defined as, 
GM | (upward) = 1/k, GM | (downward) = 1/k, 
uc а | PM . = a, : 
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A set of minimum guaranteed GM and PM may be obtained by 


defining %= min[l + g(jw)] as shown in Figure 3.3 where, 


and 





= 


Барате I2 The Nyquist Plot - Gain and Phase Margin. 
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Note that in the above definition of GM and PM, either 
$ 'or' Ki is allowed to change in the loop. The allowed 
changes are therefore very restrictive. A more useful defi- 
nition of the gain and phase margin for MIMO system that 
accounts for simultaneous changes in both Ф. and/or K, (the 
so-called universal gain and phase margin) has been derived 
moi Ret. 19]. From Figure 3.3 it can be seen that ag may 
provide a rather conservative estimate of the actual gain 


and phase margin. 





Figure 3.3 The Minimum Guaranteed GM and PM. 
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In MIMO system, gain and phase margin characterize the 
ability of the system to tolerate gain and/or phase changes 
within all loops simultaneously. Figure 3.4 shows a 


perturbed MIMO system with the assumption that L(jw) = Dia(l, 


(jw), law ) a As in the SISO case, the system will 
remain stable as longer as 1; (jw) satisfy GM INKL GM Í 
(assuming that $.s are zero ). For the case when the 
magnitude of 1(јм) are constant, the system will remain 


stable for all ?.s that 5 ۶-۷ 





Figure 3.4 The Perturbed MIMO System. 
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The set of minimum guaranteed multivariable gain and 
phase margins are determined by the minimum singular value 
(oc ) of the return difference matrix [I * G(s)], where G(s) 
is the open loop transfer matrix and consists of the plant 
and its controller. Note the similarity between the SISO 
and MIMO cases, 


sm د525‎ 141 (lL (бон) to the origin. 


о: nearness of matrix [I + G(s)] to singularity. 


Two important results that are related to multivariable 
phase and gain margin developed in [Ref. 15] are next 


presented as a theorem. 


THEOREM: The multiplicative perturbed system (Figure 
3.4) is stable if either of the following conditions hold: 


| 
۹ 


5٠۰٠۰ 2۴ ری‎ 9 (eqn 3.1) 


N 
ја 


BENNETT » ç [L (s) = T] (eqn 3.2) 


where v, c denote the maximum and minimum singular values - 
of [I*G(s)] respectively. 


It will be shown that condition 1 can be related to the 
optimality condition of LQ system and hence provides a 


useful relationship between the weighting matrices and 
robustness. 
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В. ROBUSTNESS 71 


Robustness properties pertaining to LQ system are 
closely related to the frequency domain optimality condi- 
Crone For the SISO case, Kalman [Ref. 1] showed that the 


return difference transfer function satisfies the inequality 


(Los gw) دل(‎ Ел ме (еап 3.3) 


Inspection of the Nyquist diagram in Figure 3.5 ص٦٣‎ 
indicates that а 5150 LQ state feedback has a guaranteed 
infinite upward gain margin , 0.5 downward gain margin and a 


minimum phase margin of + 60 deg. 





Figure 3.5 SISO ТО Nyquist. 
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In a similar manner, multivariable stability margin can 
be derived from the multivariable version of condition 1 in 


equation 3.1 ; 
[I + G(s)]" R [I + G(s)] > R | (eqn 3.4) 


where 
Seul -1 
G(s) - R'! B'P[sI - A] B 
апа Р > 0 satisfies the steady state Riccati equation. 
It can be shown that condition l in equation 3.1 above 
can be written as 


1/2 


go {I + R G(jw) ۷٧-٧ W1 (eqn 3.5) 


which can be reduced to 


[I + G(jwu)] > 1 (eqn 3.6) 
mivere R is diagonal ( i.e. R = pI for some positive 
Sealar °). Ke The SISO case, the multivariable LQ regu- 


lator with loop transfer matrix G(s) that satisfies equation 
3.5 and 3.6 has, in all the feedback loops, a guaranteed 


minimum gain and phase margin given by, 
GM = 1/2, and PM = +60 deg 


The case when R is not diagonal is interesting as 
condition 3.6 no longer apply and additional trade-off can 
be obtained by including off diagonal terms in the control 
input weighting matrix. This is especially useful when the 
nature and structure of the disturbances are known. Design 
involving the off diagonal R matrix will be shown in 


Appendix A. 
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IV. DESIGN PROCEDURE 


A. GENERAL 


The background material described in Chapter 2 and 3 are 
the basis for developing a computer aided design package and 
the corresponding deSign procedure. The design philosophy 
presented in this Chapter is most useful when a deSigner is 
able to characterize the desired system in terms of closed- 
loop eigenvalues and time response. Some initial physical 
insight of the weighting matrices ,their asymptotic proper- 
ties and the nature of the perturbation will be useful in 
the design process. This Chapter begins with a discussion 
of the various approaches to the LQ pole placement problem. 
The selected approach and the corresponding computer aided 
design package is presented. The design philosophy and 
design procedure for both the reduced order and full order 


model are then given. 


B. APPROACHES 


All the pole placement algorithms for LQ system that 
have been developed so far require an expression that relate 
the characteristic equation of the optimal system with the 
elements of the weighting matrices. Two such formulations 
for the stabilizeable and detectable time-invariant linear 


system (equation 2.1) and the quadratic criteria (edd 


2.3) are given by 


&(s)é(-s) = #(s)#(-s)det[I + к H(-s) Q H(s)] (en sS 


and 


D 


$ (s)$(-s) TT a eR uC QRZ (eqn 4.2) 


where ¢ (s) = det[sI -A +BF] and ¢(s) = det[sI - A] are 
the closed-loop and open-loop characteristic polynomials 
respectively. BCS = SIA) P is the open loop transfer 
matrix of the system. 

Both formulations have been used in root-locus type of 
procedure to investigate how the closed-loop poles move as 
weighting matrices changes [Refs. 5,6]. For MIMO case, 
there has been little success due to problems involving 
polynomial matrices. In this work, a different approach is 
adopted. Equation 4.l or 4.2 is formulated as a numerical 
optimization problem in which the objective function is made 


equal to the determinant part of equation 4.1 or 4.2 


| - 
Obj = det[I + RH (s) Q H(s)] (eqn 4.3) 
or 
az ٨ 
۰۱۱۰۰١۲٣٦٠٢ H (-s) Q H(s)R `] (eqn 4.4) 
For a given desired closed loop pole s = S4 , v equation 


4.3 or 4.4 becomes 
ИЕЫ РЕС заа) = о (ечп 4.5) 


Providing that Ф(5 4 )% (-sq ) is not equal to zero, the 
۲۰٠٠۲1٤ function must equal to zero if the particular Q, R 
ES to correspond to the desired closed-loop poles. 
Convergence to zero for a given set of Q and R is therefor 
automatically guaranteed. 

The pole placement problem can therefore be solved as an 
unconstrained multivariable optimization in which the 


elements of Q and R are varied to make the objective 
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function in equation 4.3 and 4.4 approach zero. This was 
done during the early phase of the work. It was later 
discovered that more insight to the problem can be obtained 
by first transforming the problem to an appropriate coordi- 
nate system and then to perform pole placement one at a time 

This is of advantage as a system designer is often satis- 
fied with several open-loop poles in a large system. 
Reassigning poles in the the reduced-order model will reduce 
computer time and memory requirement. 

The optimization routine selected for this work is the 
so-called SUMT method (Sequential Unconstrainted 
Minimization Techniques) obtained from the ADS package in 
thet sie Oe In this method, the objective function 
(eg. equation 4.3 or 4.4 ) and any constraint equation are 
formulated into an augmented objective function in which the 


problem is solved as an unconstraint optimization task. 


C: ~POLE, PLACEMENT ALGORITHMS 


For ease of implementation and better insight, only the 
problem of determining the state weighting matrix Q (given 
R) that gives a set of closed-loop eigenvalues is consid- 
ered. It must be emphasised that the algorithm can also be 
formulated to determine R (for a given Q) or to vary Q and R 
at the same time. In most cases, the present formulation is 
adequate as designers usually have some knowledge about the 
control weighting matrix. System matrix that has real and 
distinct eigenvalues is presented first, follows by cases 


where A has complex eigenvalues and repeated eigenvalues. 


1. System with Real and Distinct Eigenvalues 


— ———————M— — 


The original system given by equation 2.1 and 2.2 is 
first transformed into a diagonal form using the 


transformation given by; 
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SC) = MEZE) (eqn 4.6) 


where M is an eigenvector matrix corresponding to the system 
matrix A. and z(t) is the new state vector. The trans- 


formed system in the new coordinate is given by, 

۰ —! 

mee) = A z(t) + M Bu(t) (eqn 4.7) 
where ^ is a diagonal matrix dia[3,, 4; , 45 , ... ] 


The performance index (equation 2.3), when expressed 


in terms of the new state vector z(t) becomes, 


Ca 
I 


M х'Ох + u'Ru J dt 


К ٢ “Z u Ru) dt 


Е Í > Т" 
Л ( z Qz + u R u) dt (eqn 4.8) 


۸ 


where Q = M QM. 


It can be shown that to move an open-loop’ pole to 


^ 


its new location given by s;= sy , only Q; is required and 
other Qs ' have no effect on the pole assignment. As an 
example, to move the open-loop pole at s EX for 
NIN 245, 4, . -.-.] to its new location s*iài, , only 
О= dia[0, ES. is required. 1 


Q can then be selected according to equation 4.3 or 
4.4 using the optimization routine. As currently imple- 
mented in the program, there is no constraint equation 
formulation. Once the value of Q that satisfies the desired 
pole location is obtained, the system is transformed back to 


the original coordinate system via, 


-TA -| 


(eqn 4.9)‏ 1ھ 
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With Q known, and R given, the optimal feedback gain F can 
be obtained by solving the steady state Riccati equation as 
given in equation 2.5 and 2.6 . Since the pole placement is 
done in the decoupled coordinate system, only the eigenva- 
lues that correspond to Q.is reassigned; all other eigenva- 
lues remain unchange. It can be also shown that the 
eigenvector corresponding to an eigenvalue is also unchange, 
this property will be shown to be useful in the reduced 
order formulation of the linear quadratic problem. 

If desired, the problem here can also be formulated 
to move more than one eigenvalues in one run. This can be 
done by modifying the objective function to include more 


terms as follows; 


Obj = 


u 449 


det[I + R rs ] (eqn 4.10) 


- 


where n is the number of poles to be reassigned. 

The augmented matrix Aug = [A + BF] is NL 
computed. If desired, the procedure may be repeated to move 
other open loop eigenvalue to its specified position using 
the new Agug as the starting plant matrix. This will 
turn result in another set ot OU and d MM The effective ٣٦ ظ‎ 


and F, after n reassignments are given by 


Ө DU UA е. "ааа Im. (еча 4.11) 
апа 
Fese SEQ S She а ао ЈАВИ (eqn 4.12) 


The above pole assignment procedure can also be applied to 
an optimal system where an initial starting and R are 
given. A good example is when the control system designer 
has some knowledge of the weighting matrix but would also 


like to meet a specific time response requirement. 
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2. System with Complex Eigenvalues 


If the same similarity transformation mentioned in 
the last section is used, the transformation matrix will be 
complex. To be able to work with real matrix, an auxiliary 


transformation of the form given by equation 4.13 is used; 
x(t) = Tz(t) (eqn 4.13) 


T = ML and M is the eigenvector matrix (equation 4.6) The 


matrix L is given by, 


0١ د‎ ол ООо. 0ے‎ 
05 05] 70.0 
L = 0.0 0 1.0 
l O. 
DEO NE oe... О 
| 
(eqn 4.14) 
The transformed system is then given by, 
; -i 
moc) = A Z(t) + T Bu(t) (eqn 4.15) 
with performance index , 
J = f eTo Tz + u Ru ЕГ (eqn 4.16) 


۸ ۸ 
where Q is now given by Q = T 


0 17 
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It can be shown that to move a pair of complex 


eigenvalues given by s = a + bj, a weighting matrix of the 
form 
0 0 О 0 
^ 
0 q. 0 
^ \ ^ 
Q т 0 0 Е 
0 
0 0 
(еча 4.17) 
A A 
with q: “ Чи 1s required. 


^ 


In a similar manner, Q can be obtained by using the 
Optimization routine, with the condition q; - formulated 
as a constrained equation. Inverse transformation and 
determination of Q, and FE are identical to the 117 ول‎ 


eigenvalue case with M in equation 4.9 replaced by T. 
3. System With Repeated Eigenvalues 


In this case, the system matrix cannot be diagonal- 
ized but the general procedure given above still apply. The 
System is first transformed into the Jordan canonical Forum 


J = U À U (eqn 4.18) 


where U is a transformation matrix which is not the ip m 


vector matrix M. An example of the Jordan form is given 
below, 
À 1 0 _ 
J = |0 1 
0 0 À 
(eqn 4.19) 
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The only difference with the two procedures mentioned above 
is that the pole reassignment has to begin at the bottom of 
each Jordan block. For example, in the system given above 
(equation 4.19), the first re-assignment will result ina 


new system given by equation 4.20 


À 1 0 
ug |0 À 1 
0 0 Sd 


(eqn 4.20) 


Determination of Q and F using the optimization and the 
inverse transformation routines are identical to the 


distinct eigenvalues case with M in equation 4.9 replaced by 
U. 


Pee COMPUTER PROGRAM DESCRIPTION 


The computer aided design package developed in this 
thesis is illustrated in Figure 4.1 The program is devel- 
oped using the top-down approach with special purpose 
Subroutines called by the main control program. The driver 
program supports three independent modes with the data entry 
portion common to all modes: 

WS Data Entry : The system matrices A, B, C and/or F, Q, 
R etc are entered tnrough a qata file. Design vari- 
ables such as desired poles locations, elements of 
matrices to be varied etc, are also specified through 
the input data file. 

2. Pole Placement Mode : In this mode, an arbitrary set 
of closed-loop eigenvalues is assigned by selecting 
ehe appropriate State weighting matrix. As shown in 


Figure 4.1, the transformation matrices for various 


E) 


cases are computed first. The pole placement is 
achieved using the numerical optimization routine 
described in the last Section; Q is obtained and then 
inverse-transformed to the original co-ordinate 
system. If desired, the results can be used as an 
input to the Linear Quadratic Control Proc. ane 

3. Linear Quadratic Control Mode : This part s 
program is adopted from the OPTSYS program. Given a 
set of weighting matrices Q and R and the system 
matrices A, В апа C, it computes’ the steady state 
feedback gain F, closed-loop eigenvalues, etc. 

4. Singular Value Analysis Mode : This Portion of the 
program allows the designer to analyze various 
designs obtained from the two modes mentioned above 
in term of singular value vs frequency plot. The 
main part of this program is adopted from [Retf- 2208 

The three modes of operation mentioned above may be used 
in any order to implement specific design objectives. A 
typical design process will involve runs alternating between 
the three modes until a compromise between primary and 
secondary design objectives is achieved. Record of a 
typical design run together with a complete listing of the 
main program and their non-standard subroutine are given in 


Appendix B and C. 


E. DESIGN ۳ 


The Linear Quadratic constant state feedback design 
philosophy for the linear time-invariant model is illus- 
trated in Figure 4.2. It assumes that the location of the 


eigenvalues and system time response are the main design 


objectives. These objectives can be achieved using the pole 
placement procedure developed in this thesis. Single or 
multiple reassignment can be made іп one run. A number of 
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Computer Aided Design Package Organization. 
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designs can then be obtained using different starting 


control weighting matrices, different state weighting 
matrices and different assignment sequences. Physical 
constraints such as control input amplitude, control ٦ 


energy as well as general system properties such as asymp- 
totic behavior are heavily relied upon during this process. 
After the major objectives are satisfactorily achieved, 
secondary design objectives are considered. These include 
feedback gain reduction (by increasing the control weighting 
matrix), robustness (in terms of minimization of system 
Sensitivity to modelling errors and/or parameter varia- 
tions), zeros locations, eigenvectors assignments. The 
extra degrees of freedom available in the MIMO state feed- 
back system can often provide a means to improve these 
secondary objectives while only slightly modifying the 


closed-loop pole assignment and thus the time response. 
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V. DESIGN 8 5 


The design procedure described in Chapter 4 is illus- 
trated in this Chapter by two design problems. A typical LQ 
structure and its properties, together with the pole place- 
ment procedure is first presented with a 2x2 model. A prac- 
tical design problem is then presented for the highly 
coupled lateral channels of a CH-47 Helicopter. The 
resulting LQ designs are compared with other multivariable 
state feedback designs [Refs. 2232208 It is shown tma 
the procedure developed here is a viable tool for robust 


constant feedback controller design. 


A. INTRODUCTORY 2X2 PROBEEM 


This problem formulated in reference 23 serves to demon- 
strate how a highly cross-coupled multivariable control 
problem can Бе formulated and solved as a linear quadratic 
design problem, using the pole placement procedure. The 
problem provides excellent insight into the structure of the 
multivariable LQ system and its built-in robustness Со 
cross-coupled perturbation. 

Figure 5.1 shows a diagram of this basic 2X2 system in 


which the plant is given by 


x(t) = Ax(t) + Bu(t) (eqn 5.1) 
ОЕ TCE) (eqn 5.2) 
where 
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Figure 5.1 





A Simple 2X2 MIMO Model. 
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This system has open-loop eigenvalue at -l, -1 and i 
therefore stable. The b, term in the control matrix B is 
purposely made large to produce the cross-coupled effect 
from channel two to channel one( see Figure 5.1). The 
design requirement is to select a set of feedback gain F 
such that the closed-loop eigenvalues are at  -2, -2. 
Assuming that this is the only requirement, it 15 not ٢ لګ‎ 


cult to see that a unity feedback law of the form (Figure 
5 295 


1 | ۱ 


u, - 1 0 X, ll, 
z F | 
u, 0 7 X. Ча | 


(eqn 5.3) 


will produce the desired closed-loop system given by, 


x 0 SOE 1 50 u | 


(eqn 5749) 


The above design seems to be acceptable as far as eigen- 
values or time response is concerned. It is now shown that 


when robustness of the system is considered, the unity feed- 


back gain controller performs rather poorly. On the GENES 
hand, design using LQ Pole Placement type of formulation 
will result in robust controllers: To demonstrate the lack 
of robustness of the unity feedback design, the feedback 


gain matrix F as given in equation 5.4 is perturbed slightly 


(by +5%) and the eigenvalues of the resulting closed los 
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Figure 5.2 Unity Feedback for the 2X2 Model. 


System matrix are calculated. The absolute and percentage 
errors inthe eigenvalue due to the perturbation of the 


closed-loop system are given in Table I 


It can be seen from the above that the first cut design 
is very susceptible to model and feedback perturbation. A 
5 % changes in the f,, term has resulted in a large shift 
( 1605 ) in the closed-loop eigenvalue. Lack of robustness 
in the unity feedback design can also be seen in terms of 


the minimum singular value plot of the return difference 


matrix in the frequency domain. For the unity feedback gain 
system shown in Figure 5.2, the return difference matrix is 
given by, 
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TABLE I 
PERTURBED EIGENVALUE FOR 2X2 MODBEE 


Perturbation Absolute Changes Percent Changes 
ID RS 52) in Eigenvalues in Eigenvalues 
f, and faa | 0.05 ٣۶۶ص‎ 2.5 2 ш 
122 0.0 TT OTO 0.0 OD 

Ёз | 0.67. ؤؤ ؤ ؤ ؤ‎ 11.5 ICD 

А11 رط‎ 1658+ 37. ٧ ؛‎ ٣۳ 
ЗЛИ) SO san 
I + G(s) = 
0 Е) 
E 
(eqn 5.5) 
where G(s) eee ean is the loop transfer matrix as ае 


cated in Figure 5.3 

The multivariable Nyquist diagram ( locus of det[I+G(s)] 
for the system is shown in Figure 5.4 If this diagram is 
interpreted as for a single input system, the (-1/2,9) gain 
margin and (+160) phase margin would lead one to believe 
that the design is a good one . This has been shown to be 
not the case, a 5% perturbation in F would cause the system 
to become unstable. The above clearly demonstrate the 
inadequacy of the classical method in evaluating stability 
margin for MIMO system. 

Robustness properties of the unity feedback controller 


Will now be analyzed in terms of singular value as discussed 
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Figure 5.3 boop ttanster Matrix = F = L. 


in Chapter 3. The minimum singular value plot of [I+G(s)], 
BEN Rhown in Figure 5.5 as a function of frequency. The lack 
of robustness is clearly indicated by the relative small 
singular value at frequency around 2 rad/s. Using the 
universal phase and gain margin chart developed in 
meet. 19], the minimum singular value at this frequency 
etam 25 5691945 to a gain margin of (0.91, 1.0) and a phase 
margin of (+ 4 deg). 

IS now shown that formulation using LQ approach and 
the pole placement procedure developed here will result in 
robust design that meet the time response requirement. 
Furthermore, better insight of the design process can be 
obtained from the procedure to be described here. The first 


step in the LQ design is to determine the asymptotic 
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Figure 5.4 Multivariable Nyquasterlor 


properties of the system, i.e. movement of closed-loop poles 
as R 40. Using results from Chapter 2, it can be estab- 
lished that as R increase from O toe, both closed-loop 
poles move from infinity on the real axis to the open-loop 
poles location. None of the closed-loop poles stay finite 
as R»O since the dimension of the input control vector is 
equal to the dimension of the state. Assuming that R = 1 
the pole placement is accomplished in two steps. The fir < 
step is to move the open-loop pole at -l to -2.0. As the 
system matrix À given is already in Jordan form, no trams. 
formation is required. The pole placement program puts the 
pole at -1.9987 With, 
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0.001 0 0 0 


0 0 0.01983 0 


and 
1.9915 0 
А суз 
0.01983 10 
In the second step, the other open-loop pole at -l. is 
moved to -2.0 , using Agyg as the new plant matrix. The 
resulting Q and F and the augmented plant matrix become, 
= 7 -149.9 0.9994 249. 970 
جک‎ 5 
-149.9 74993 -0.00841 0.4211 
Ine eftective O und Fe required to move both open- 
1оор ро1еѕ ас -1.0 ёо -2. аге О = 0 + О, апан 0 
Ез аз shown in equation 5.6 bellow, The pole placement 


procedure is completed with the final eigenvalues placed at 
17آ 1 91 ٴ ٴ"۰)‎ 7۶٥ 


200915; د یں‎ 0.9994 -49.978 
Q, 7 Е, = 
-1۹49 9 250905) 0.0114 0.4211 
ES L. 
(eqn 5.6) 


p 


The singular value plots of the case where R = I 
together with the unity feedback gain ( non-LQ design ) are 
shown in Figure 5.6 The well established fact that LQ 
design possesses ( 1/2,9 ) gain margin and (+#60deg) phase 
margin can also be readily observed from the same figure as 
the minimum singular values of [I*G], “ , is greater than 
one for all frequency. Changes in closed-loop eigenvalue 
for a small perturbation in F is again computed as shown in 
Table II . It can be seen that the LQ design is robust with 
the largest percentage change in eigenvalues location of 
only 10%, when compared with the 160% change in the unity 
feedback design. 


TABLE II 
PERTURBED EIGENVALUE FOR LQ DESIGN 


-о æ ма мом ее ес е мое е см е е ее ао еее а е е е се е = ч = е еко ее се е се се е се е ем پس‎ е ې‎ оо в аә بس‎ е е ې‎ мм 


Perturbation Absolute Changes Percent Changes 
meer (+54) in Eigenvalues in Eigenvalues 
ENSE. | 0.2149 , 0.142 | 10.7 ,7. ` 
ШЕ са ОД. DNO OS Mh, ЕРКО 
е ү (00592, 00578 Жү ор, 29 , 
جک‎ | ПО, ооо" у DONE ۲ 


=» وس‎ е е см ее پس‎ ме ом سا‎ ес мо ме е ee оо мо ес мм ее сме оо мо мо то е е پا‎ мо ہے پا جسہ سیا ېی سی پس‎ мо ج‎ мо е ہے کے چوس‎ е смо оо س جس‎ 


The pole-zero plots of various closed-loop transfer 
functions of the closed-loop transfer matrix for both the 
unity feedback and LQ design are compared in Figure 5.7 and 
DS . For the unity feedback design, zero at minus three 


for the input 2 to output 1 channel corresponds to the 
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minimum singular value frequency ( point À in Figure 5.6). 
It is in fact this zero that causes the the system to be 
sensitive to perturbation. For the LQ design, the built in 
robustness cause the zero at minus three to move toward the 
pole location at -2. 

Improvement in the robustness from LQ design is now 
analyzed in terms of Bode plots. The open-loop gain and 


phase vs frequency plots for a MIMO system can be obtained 


wom the open-loop transfer function matrix, G(s). For full 
State feedback, de scc e (if the open-loop plots 
without feedback is considered, G(s)=C(sI-A) İB ) The 
matrix G is a matrix-valued rational function of s, it 


describes how the system (with or w/o feedback) appears to 
1ts environment. It is an external description of the 
system and is closely related to the zeros of the system. 
The open-loop Bode plots for the two designs are compared in 
Figures 5.9 to 5.11 As b2, =0, there is no coupling from 
channel 1 to channel 2. All channels have a -20 dB/decade 
slope at high frequency which is in agreement with earlier 
observation that each channels has опе finite zero. It is 
interesting to note that the unity feedback and LQ design 


result in almost identical gain vs frequency plots for 


direct channels (i.e. channel 1-1 and 2-2). Any classical 
Single loop type of analysis will not be able to detect any 
difference between the two designs. Оп the other hand, 


cross-coupling effect can be readily seen from the gain vs 
Breguency plot in channel 2-1 (Figure 5.10). The unity 
feedback design is characterized by the rather large channel 
lI Bain at low frequency. A very small perturbation in 
channel one's parameter can change the system behavior 
considerably. This has been illustrated earlier Бу 
perturbing the feedback gain. It can be seen from the 


figure that LQ design reduces the crossfeed gain by large 
amount. 
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2۰۰۰۰۸۶ HELICOPTER DESIGN PROBLEM 


The design procedure described in Chapter 4 is further 
illustrated in this section by an actual design problem. A 
controller is designed using the Linear EC Pole 
Placement approach for the lateral dynamic model of a CH-47 
helicopter. The resulting design is then compared to multi- 
variable state feedback controller given in [Ref. 23]. 

The highly coupled two inputs lateral axis model of the 
CH-47 helicopter is used as a full order system. The state 


vector x(t) and control input vector u(t) are given by, 


xX, =v = Y-axis velocity (ft/sec) 

X2 = p = Roll rate (rad/sec) 

X, = r= Yaw rate (rad/sec) 

x, = $ = Roll angle (гаа) 

u, = 8y, = Yaw rate rotor deflection control (inches. ) 


uy е“ дес = Roll rate rotor deflection control (inches. ) 


The state variables and the body axes of the aircraft are 


illustrated in Figure 5.12 The yaw and roll rotor deflec- 
tion control produce changes in the yaw rate, slide slip 
angle, roll rate and bank angle. Assuming full state feed- 


back, the A, B, and C system matrices are given by, 


E227 1.420 SUNS Dog 


À = . 01 =0.7 307 07 0.0 
0.04 - 0.05 7-05 oF 
o: L- 0.11 0. 
0.12 095 1 0 0 0 
> 0.04 807 Cm = - 0 1 О О 
. 34 .020 0 0 1 0 
0.0 ‚0 | 0 0 0 1 
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The open loop eigenvalues of this system are; 


OF 40.2065 
5, = -0.0503 
S 3 = -1.0498 
pee = -2.1263 


The open-loop system is not stable and the time response 
is shown in Figure 5.13 for zero input and an initial condi- 
tion of Ф(0 )=0.1 гаа 

The design requirements are to satisfy specification in 
terms of the step input response of the roll attitude 
channel CPE J: Stability margin requirements are those 
given by standard military specification. In [Ref. 23] 
three designs were obtained to satisfy the desired time 
response performance specification. All three control laws 


are of the form given by, 
u(t) = -Fx(t) + M) (eqn 5.7) 


The values of F and h are summaried in Table III . It was 
shown in [Ref. 23] that two of the designs (design 1 and 2 ) 
were extremely sensitive to model errors and perturbations. 
It is now shown that LQ formulation using the pole placement 
procedure developed here will result in robustness design 
and yet satisfy the conventional time response criteria. 

It is assumed that the closed-loop poles requirement are 
the same as those obtained in the AlphaTech's design |. 
۳-12, -12.51, -9.652, -2.125). The first step in the 
design is to establish some asymptotic properties of the 
System. The dimension of the control input (2) is less than 
the dimension of the state(4). As shown in Chapter 2, this 
implies that at least two closed-loop poles approach minus 
Eunity in the complex plane when Е ә 0. When R-~ 0 


(meaning no control input allowed), the closed-loop poles 
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TABLE III 
ALPHATECH DESIGNS ( 1, 2, AND 3) 


Designs Feedback Gains (F) h 
Mone | -1.72 -23.5 70.6 595 | 595 
0 ЕО 0) 168 ۶ 9 27.99 
two 1 5 18.3 142,0 142-0 
“0. 01 -1.592 -0.189 -1.47 = А] 
three 0 0 55 0 
0 -4 0 - 27 1 =] 


approach the open-loop poles or their mirror images if they 
are in the right hand plane. e. -0.2065, -0.0503, 
04987, and -2.12635). 

The next step in the design procedure is to select a 
Suitable starting control weighting matrix. Based on guide- 
Mtesgiven in [Ref. 4] and assuming that u,,a, and u,,,, are 
equal to l inch. R = [I is selected for the initial design. 
(In the present formulation p is not needed as only Q is 


varied to place the pole ). 


The sequence of the reassignment is determined next. At 
present there is no known established guideline for 
selecting the prefered sequence of moving the poles. Two 


extreme sequences are considered as follows, 
A. Move the left most open-loop pole to the left most 
desired closed-loop pole etc. The reassignment 


sequence then become: 


I move pole at 2.1265 to -25.12 
1076 pole at 1.048) to -12.51 
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3. move pole at -0.2065 tom. n 
4, move pole at -0.0503 to -2.125 


B. Move the right most open-loop pole to the left most 
desired closed-loop pole etc. For the problem given, 


the reassignment Sequence are 


1 move pole at +05 2055 11 7 
2 move pole at -0.0503 to -12.51 
3. move pole at 10490 +76 
4 move pole at -2.12 TEGO О 
(close, so no move required) 

Using the pole placement program developed in this 
thesis, the corresponding Q and F for each of the reassign- 
ments are obtained and tabulated in Tables IV and V for the 
two cases selected (design LQ-A and  LQ-B). It can be seen 
from the table that difference reassignment sequence results 
in different set on and pP E In general, once the 
matrix R and the n closed-loop poles are selected, there are 
n(n*1l)/2 extra degrees of freedom available in Q. This als 
verifies the well-known fact that the general solutions of Q 
and R for a given set of closed-loop eigenvalue are non- 
unique. The elements of the Q matrix obtained during each 
reassignment depends on the transformation matrix (M, L or 


U) and hence the eigenvectors that are uSed to construct 


them. The important of eigenvector type of assignment is 
evident as Q, F and the resulting closed-loop poles depend 
on the eigenvectors used. The designer can shape the 


design by choosing appropriate eigenvectors for the trans- 
formation matrix. The application of these extra degrees of 
freedom will be discussed in the following section. The 
resulting designs (LQ-A and LQ-B) are now compared with the 
first desienm in 6 299 
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TABLE IV 
RESULTS TRON POLE PLACEMENT SEQUENCE(LQ-A) 


Q and F obtained during each reassignment 
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TABLE V 
RESULTS FROM POLE PLACEMENT SEOUBENGCPIENEI 
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the yaw rate response; the yaw rate response is well damped. 
The LQ-B design appears to be a better design as coupling 
among modes are small and can not be detected from the time 
response plot. It is noted that the difference in response 
for the same set of closed-loop eigenvalue 15 due to 
coupling among various modes through their respective eigen- 
vectors. The ‘overshoot’ in this case is obviously not due 
to complex conjugate poles as all closed-loop poles in the 
three designs are on the real axis. The set of final eigen- 
vector for the three designs are tabulated in Table VI 

Inter-mode coupling for the yaw rate response in the 


AlphaTech design and roll attitude response in the LQ design 


case A can be readily seen from the table. The issue of 
eigenvector assignment will be discussed in the next 
section. 

TABLE VI 


CLOSED-LOOP EIGENVECTORS 


ae 79907 See 22 299724503 -11.9449 

aipha] -0.39336 USA -0.48016 -0.41854 
one 0.107067 0.00822 -0.32958 -0.33640 
=0. 21300 0010955 O 61191 20 ig 
-0.00029 -0.00488 0.04306 00593 

! е -9.26645 - 1271520065 -25%200990 
LQ-A 0000 0 -0.15459 .01894 20.887053 
0.18018 0.008218 0.99664 0 28 

m0). 9560 3 pO 01258 0.00396 0.00031 

=02 05029 0.10614 0.07964 707002 
295.17894 295627038 -12.51644 22857630 

IFO-B S0 12652 2 +0 19209 0.05042 UE ОЛУ 
0.99118 "0 5055 10.9015 ) 07 
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Figure 5.14 
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Figure 5.15 Closed-Loop Time Response Plot (LQ-A). 
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Closed-Loop Time Response Plot (LQ-B). 


Figure 5.16 
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The difference in time response performance for the 
three designs can also be interpreted in term of the closed- 


loop pole-zero plots in the frequency domain. The pole-zero 


plots for the designs, with $c аз the input, are shown in 
٦۔٥‎ 5.17 to 5.19 . AS all pole locations are the same 
for the three designs, they are indicated only once in each 
diagram. 


In the AlphaTech's design, the zero locations for the 
various channels clearly indicate the weakness in the 
design. Most transmissions zeros are located away from the 
poles" locations. The system is therefore strongly coupled 
to its external environment. А good example is the lightly 
rupaq” pair of zero at (-1.05,+7.31j) for the v-%@ chan- 
nels. It is in fact these undesired zeros that reduce the 
overall robustness of the system. The mechanism of robust- 
hess improvement in the LQ design can also be seen in the 
pole-zero plots in Figures 5.18 and 5.19 The built-in 
robustness in the LQ design causes the zeros at various 
channels to move to locations where their transmission prop- 
erties can be canceled by the closed-loop poles. An excel- 


Eu example is-in-—the v- é&. channel, where the zeros at -2, 


Er and the mirror image of +24.0 are fairly close to the 
Suesed-loop pole locations at -2.12,  -9.26 and -25.12. In 
the case where the zero from LQ design are not close to the 
Sed loop poles-(r- e$. channel for LQ-A ) the zeros are 
"well damped" and their transmission properties can be 
neglected. 


Robustness properties of the three designs presented 
here can also be analyzed from the open-loop’ Bode plots. 
EE open-loop transfer function gain plots for channel 1-1, 
|| 2-1 and 2-2 for the three designs are shown in Figures 
END to 5.24 Cross coupling problem for the AlphaTech 


design is clearly indicated by the relatively high gain of 
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the channel from input 2 to input 1 (2-1). Both LQ designs 


presented here reduced the gain in this channel by more than 


40 dB inthe frequencies of interest. The bandwidth is 
reduced from above 100 rad/s to about 8 rad/s. Unlike the 
simple (2x2) System where there was little difference 


between the direct channels of the two designs, gain adjust- 
ment is observed in all channels. For example, LQ designs 
reduce the gain in channel 1-l but increase the gain іп 
channel 2-2. There is also a slight increase in the 
Eoupling channel 1-2. The overall effect is that of gain 
balancing, gains in channels that are affected by  cross- 
coupling perturbation are lowered together with some adjust- 
ments in other channels. Different  reassignment sequence 
results in different adjustments. The designer has to 
choose a set of gain curve depending on the particular 
requirement. The relationship between the open-loop Bode 
plots and the zeros of the system is also evident from these 
diagrams. The low frequency resonance (or ‘peak') in the 


gain vs frequency plot for the AlphaTech design correspond 


to the undesirable Zero mentioned earlier. In the LQ 
designs, these low frequency resonances are absent because 
of the more desirable zero locations. It is also inter- 


esting to note that the LQ-B gain curves fall nicely in 
between the other two designs. This fact, together with its 
Detter overall low frequency gain characteristic, may 
account for the better time response behavior of the LQ-B 


design. 


Exc final comparison, the closed-loop frequency 
response plots of the various channel for the three designs 
u shown in Figures 5.25 to 5.32 . At high frequency, all 
the gain plots approach the -20dB/decade slope as all eigen- 
values are on the real axis. The AlphaTech design has a 


near Odb gain for frequency up to about 2 rad/s. It appears 
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to be an acceptable desipn but was shown to have undesired 
transmission zeros in other channels; this can also be seen 
from the closed-loop plot for the ro channel in Figure 
5.31 . Design LQ-A, obtained from the reassignment sequence 
A, is characterized by very low DC gain (-9dB) and ‘peak in 
frequencies near 10 rad/s. These have been shown, both in 
the time response and pole-zero discussions that it will 
have undesired effect on Спе рос Ера ЕЕ Design LQ-B, 
obtained from the second reassignment sequence appears to be 
the best compromise. All channels (Figures 5.30 to لر‎ 
have flat low frequency characteristic and coupling between 


modes are almost absent. 


C. DISCUSSION AND CONCLUDING REMARKS 


A new computer aided design procedure for the multivari- 
able linear time-invariant system using Linear Quadratic 
Pole Placement formulation is presented. The two design 
examples presented above have served to demonstrate the 
complexity involved ina multivariable design. The main 
problem lies in the fact that the solution of a MIMO problem 
is in general non-unique. This was shown in the helicopter 
problem where different approaches result in different 
designs, although the closed-loop eigenvalues for all 
designs were the same. It was also shown that the ех 
degrees of freedom in MIMO system design can be accounted 
for by analyzing the singular value plots, transmission zero 
movements and closed-loop eigenvectors of the designs. It 
becomes apparent that the success of any multivariable 
design methodology hinges on how to make use of these extra 
degrees of freedom available in MIMO system. It was also 
shown in the above design examples that the Linear Quadratic 
Pole Placement procedures possess such quality. Some unique 


properties of the method presented here are as follows; 
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Figure 5.24 Closed-Loop Bode Plot- AlphaTech 1-1. 
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Figure 5.28 Closed-Loop Bode Plot- LQ-A 1-3. 
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Figure 5.29 Closed-Loop Bode Plot- LQ-A 1-4. 
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Figure 5.30 Closed-Loop Bode Plot- LQ-B 1-1. 
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Figure 5.32 Closed-Loop Bode Plot- LQ-B 1-4. 
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It guarantees stability 210 ٢ 

It has partly overcome the main difficulty لر‎ 
applying optimal control theory to practical r 
lems; і.е. that of selecting suitable performance 
indices (О and R matrices). The designer can now 
choose to compute Q (if he knows somethings about R) 
or vice versa or compute both Q and R if both of them 
are unknown. 

3. It has built-in robustness to model's error ame 
perturbation. 

4.  Reduced-order type of problem formulation is possible 
as the procedure can reassign poles опе ас a time 
without affecting others. 

On the other hand, the design procedure developed here 
is not without some shortcomings. Three main problem areas 
are described below, one of which can be overcome by addi- 
tional programming efforts while the others are still active 
research areas pursued by many researchers. 

1.  Eigenvectors Assignments: It was shown in the design 
examples presented above that while the overall speed 
of response of the closed-loop system is determined 
by its eigenvalues, the ‘shape’ of the transient 
response depends on the closed-loop eigenvectors. 


The problem of eigenvectors assignments was first 


formulated in RP SU 1! Since then many eigenvec- 
tors assignment algorithms have been developed 
[Refs. 17,24] for use in multivariable design. In 


principle, the eigenvectors assignment routine can be 
incorporated. 5 the coordinate transftformat im 
portion of the OPTPP program (as indicated in Figure 
уел aor In essence, a new Similarity 
transformation different from M in Equation ٥٣ 
computed once the eigenvalues and eigenvectors are 


Specified. As most eigenvector algorithms available 
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at present are iterative in nature and their inclu- 
Sion requires major programing efforts, it is recom- 
manded for future work. 

Perturbation or Model's Error Structure: The guaran- 
teed robustness obtained from the LQ formulation 
given by equations 3.4 and 3.5 ensures that the 
perturbation or model's error (L(jw)) is sufficiently 
small so that the closed-loop system remains stable. 
However, the above only applies to simple model's 
error structure where both L(jw) and R are diagonal 
matrices. There may be cases when the above equa- 
tions do not hold and hence the design becomes very 
conservative. An example of this nature is shown in 
[Ref. 19]. An ad hoc solution is to use non-diagonal 
control weighting matrix as mentioned in Chapter 4. 
Two designs for the helicopter problem are obtained 
using  non-diagonal R. Their results in terms of 
Басат ше plots, Bode plots etc are compared 
with other designs in Appendix A 

Reassignment Sequence: As shown in the last section, 
different reassignment sequences result in different 
designs. At present there are no known methods for 


determining the best reassignment Sequences. 


МЕЈ 


VI. CONCLUSIONS 


It was demonstrated that the general pole assignment 
problem in multivariable state feedback control system 
design can be formulated using the Linear Quadratic Control 
approach. This method of formulation is effective for two 
main reasons; First, the extra degrees of freedom available 
in a multivariable system structure is utilized to produce 
designs that are robust to perturbations in the system and 
gain matrices. Secondly, the classical difficulty ties 
selecting Suitable performance index in optimal control 
formulation was partly overcome, as designers now have the 
flexibility of specifying only Q (the state weighting 
matrix) or R (the control weighting matrix) or both QU ande 
as the design parameters to be varied. In other words 
knowledge of the performance index which ideally should come 
from physical argument is used to the best of designer's 
advantage. In addition, the structure of the present formu- 
lation is such that eigenvector assignment can be further 
incorporated into the procedure. The above properties, when 
combined with the reduced-order formulation capability, have 
been shown to be very versatile and have important impact on 
the performance of the resulting design. 

The optimal pole placement (OPTPP) program developed 
here is combined with other well established routines to 
fomm a computer aided design and synthesis package. 
Together with the design procedure and philosophy presented 
here, it provides the control system designer an excellent 
and viable tool to solve complex multivariable problems. 

The procedure was applied to practical test examples, 
and numerical results were presented and discussed. Results 


indicated that all controllers obtained from the formulate 
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given here were stable and robust. Introducing perturbation 
in the system matrices leaded only to small errors in the 
assigned poles. The main shortcoming of the design proce- 
dure 15 the ad hoc nature in which the poles are reassigned. 
More research is required to develop a systematic way of 
assigning poles and its eigenvectors thus allowing the 


designer to optimaly shape the response of the system. 


ЭЭ 


APPENDIX A 
NON-DIAGONAL R DESIGNS 


The two design examples presented in Chapter 5 were 
based on diagonal control weighting matrix )1 ۶٥دد‎ ROM It 
is now illustrated that design using  non-diagonal R will 
provide yet another degrees of freedom available to the 
designer. This type of formulation is especially useful 
when the model's or error structure is known. For multipli- 
cative type of perturbation, the effect of control weighting 
matrix on the system stability has been explored in 
[Ref. 12]. In essence, the Selection of the R matrix deter- 


mines the coordinate frame in which the sensitivity assess- 


ment is to be made. This can be readily seen from equation 
3.1, the general form in which R is non-diagonal 15 БАШ 
below; 

0 (> —1 - 

СРС ٣ = F] <+ e TI (9 (eqn A.1) 


In general, the sensitivity of the system to perturba- 
tion ina particular L(s)'s direction can bé reduced 
making = ہہ‎ ٤+٥٣ S smale This can be done 
simply by choosing an appropriate R. The main problem with 
this kind of approach is that L(s) must be known 26ت‎ ٣ 
for all frequency for the computation "of з ۶۴نھ]‎ ٠ (s R ОШ 
E In addition, the worst case direction of L(s) الک ال‎ 
known otherwise the resulting design may be too 
conservative. 

The helicopter problem presented in Chapter 5 is now 
analyzed using the non-diagonal control weighting matrix. 


Two designs are obtained as follows; 
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Design one (LQ-C): 


Design LQ-C is obtained using an assignment sequence 
similar to LQ-A with no off-diagonal element їп К. Design 
LQ-D is obtained using placement sequence for design LQ-B 
but has off-diagonal element in R. Q and F obtained during 
each move and the final Q, and Fẹ, are tabulated in Tables 
VII and VIII. Singular value plots and the open-loop Bode 
Plots are compared with other designs in Figures A.l to A.5 

Both design has singular value greater than one for all 
frequency. It is interesting to note that similar reassign- 
ment Sequence produce similar singular value plots. This 
can be readily seen by comparing LQ-A with LQ-C and LQ-B 
with LQ-D in the figures. The effect of using non-diagonal 
R merely changes the shape of the singular value plots. For 
the cases presented here, both non-diagonal R designs are 
Slightly less conservative (having lower singular value). 

As in the Singular value plots, the shape of the open- 
loop Bode plots (with feedback) are closely related to the 
reassignment Sequence.(compare Bode plots of LQ-A with LQ-C 
EXUGEO-B with LQ-D in Figures A.2 to A.5 . Results from 
Ele zero maps also indicate similar trends. 

In summary, it is demonstrated that designs using non- 
diagonal control weighting matrix (including possibly off- 


diagonal elements) Provide yet another means of ‘fine 
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tuning’ the design. This capability, such а5 1۹11 o sS 
diagonal elements directly in design, is unique to the 
present formulation. Robustness between the upper and lower 
crossfeed (of L(s)) can be controlled by adjusting the eum 
tive weighting of the upper and lower (or off-diagonal ) 
elements of the control weighting matrix. It must be empha- 
Sized that this kind of fine tuning is only possible for a 
class of rather well-defined structure of L. In practices, 
other constraints, such as the energy of the control input, 
the conditioning of R etc must be considered over the range 
of frequencies. The key issue of model's error structure 
and how it can be used in multivariable control system 


design iS currently being investigated by many researchers. 
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TABLE VII 
RESULTS FROM POLE PLACEMENT SEQUENCE (LQ-C) 
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